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We examine the step dynamics in a 1+1 dimensional model of 
epitaxial growth based on the BCF-theory. The model takes an- 
alytically into account the diffusion of adatoms, an incorporation 
mechanism and an Ehrlich-Schwoebel barrier at step edges. We 
find that the formation of mounds with a stable slope is closely re- 
lated to the presence of an incorporation mechanism. We confirm 
this finding using a Solid-On-Solid model in 2+1 dimensions. In 
the case of an infinite step edge barrier we are able to calculate 
the saturation profile analytically. Without incorporation but 
with inclusion of desorption and detachment we find a critical 
flux for instable growth but no slope selection. In particular, we 
show that the temperature dependence of the selected slope is 
solely determined by the Ehrlich-Schwoebel barrier which opens 
a new possibility in order to measure this fundamental barrier in 
experiments. 

PACS: 81.10.AjTheory and models of crystal growth; physics of crystal 
growth, crystal morphology and orientation 

1 Introduction 

Molecular beam epitaxy (MBE) has attracted much interest from both, theo- 
retical and experimental physicists. On the one hand it allows the fabrication 
of high quality crystals with arbitrary composition and modulated structures 
with atomically controlled thickness ][]]. On the other hand it represents a 
model of nonequilibrium physics which still lacks a general theory [fj] . In par- 
ticular, the appearance and the dynamics of three dimensional (3D) struc- 
tures (pyramids or mounds) in crystal growth are not well understood in 
terms of the underlying microscopic processes. 
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A long time ago Burton, Cabrera and Frank introduced the BCF-theory of 
crystal growth ||. Within this theoretical approach the crystal surface is 
described by steps of single monolayer height. The evolution of the surface 
is calculated by solving the diffusion equation on each terrace. Within this 
framework the growth of spirals and the step flow has been investigated. Elk- 
inani and Villain investigated such a model including the nucleation probabil- 
ity of new islands @] . They found that the resulting structures are unstable. 
Towers appear which keep their lateral extension and grow in height only. 
They called this effect the Zeno-effect. The same observation has been made 
with a "minimal model" of MBE where fast diffusion together with a high 
Ehrlich-Schwoebel barrier has been implemented ||. 

Even though the Zeno-effect has been observed recently on Pt(lll) ||, quite 
typically a coarsening process with appearance of slope selection emerges 
which has been reported for such diverse systems as Fe(OOl) |7], ||, Cu(OOl) 
@, [Kj, GaAs(OOl) 0, and HgTe(OOl) jLj. In addition, slope selection 



seems to be the generic case of solid-on-solid computer simulations [14], [12 
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In terms of continuum equations the selection of a stable slope has been 
related to the compensation of uphill and downhill currents |16|, [XT]. An 
uphill current can be generated by an Ehrlich-Schwoebel barrier |18|, |19[ . The 
barrier hinders adatoms to jump down a step edge. Hence, more particles 
attach to the upper step edge which leads to a growth-instability PU| and 
3D-growth. 

Another process, which constitutes a downhill current, has been recognized 
using molecular dynamics simulations [^T|, [2^]. Such diverse mechanisms as 
downward funneling, transient diffusion or a knockout process at step edges 
lead to the incorporation of arriving particles at the lower side of the step 
edge. In addition, it has been suggested that such a process is responsible 
for reentrant layer-by-layer growth |p3| 



Recently we have proposed a simplified model of epitaxial growth quite sim- 
ilar to the "minimal model" of Krug p4| . In particular we found that an 
incorporation mechanism is crucial to achieve slope selection. However, one 
simplifying assumption of the model is an infinite Ehrlich-Schwoebel barrier. 

In this article we will present in more detail the argument leading to slope 
selection and we will generalize our results using a continuous step dynamics 
model analogous to |0J. In sec. |2] we will introduce our extension of the BCF 
theory and will discuss the relation to existing results (sec. || and[3|). Typical 
mound morphologies and the growth dynamics are compared in section 01. 
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Figure 1: Density of diffusing adatoms on a terrace of size I. The origin of 
the x-axis is chosen to be in the middle of the terrace. 

Afterwards we will investigate the emergence of slope selection within the 
framework of this model (sec. |5|). We will show that the selected slope 
has a temperature-dependence which is solely determined by the Ehrlich- 
Schwoebel barrier. Hence, the determination of selected terrace widths in 
experiments would give direct insight into microscopic properties such as 
the Ehrlich-Schwoebel barrier. We confirm the predicted importance of the 
incorporation mechanism using a kinetic Monte-Carlo simulation of a Solid- 
On-Solid model in sec. |6|. Another effective downward current could be due 
to detachment from steps and subsequent desorption. We will show in sec. [7| 
that slope selection cannot be achieved by these two processes alone. In 
section ^| we will calculate the saturation profile in the limiting case of an 
infinite Ehrlich-Schwoebel barrier. 



2 BCF theory 

The model is based on the Burton-Cabrera-Frank model in 1+1 dimensions. 
Within this framework the crystal surface is specified by the position and 
direction (upward or downward) of steps. Figure [I] shows the crystal surface 
from the point of view of the BCF-theory. It is a coarse grained view - the 
detailed positions of atoms are not important. However, the terraces of the 
height of one atomic monolayer (ML) can still be distinguished. The most 
fundamental assumption is that at each time t the adatom concentration 
p is a function of the step positions only. In other words, the diffusion of 
adatoms is considerably faster than the step velocity. Thus, the diffusion 
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equation becomes 

^(x,t) = = DV 2 p(x,t) + ^ (1) 

where D is the diffusion constant and F j a is the flux density with a denoting 
the lattice constant. Hence, 1/F is the time necessary in order to deposit 
one monolayer. Up to now, this equation was solved with special boundary 
conditions at x = —1/2 and +£/2 in the literature. Clearly, the boundary 
conditions are chosen depending on whether the terrace is a vicinal, a top, 
or a bottom terrace. In the following we will discuss the typical case of a 
vicinal terrace. The extension to top and bottom terraces is straightforward. 

To include an incorporation mechanism it is necessary to extend the theory. 
We assume that there exists an incorporation radius such that all particles 
arriving close to a downward step within this radius immediately jump down 
the step edge. Hence, one has to split the density of diffusing particles into 
two regions. The first region close to the upper edge where eq. (|I|) holds, and 
the second one given by the incorporation radius close to the downward step 
where no particles arrive (F = 0). To describe the motion of steps the flux 
of incorporated particles must be taken into account separately. 

In the following we will discuss in detail the situation £ > R inc as sketched 
in fig. p]. For smaller terraces only one region exist and the calculations are 
much easier. Since our analytical calculations will show that £ > R\ nc is the 
generic case we concentrate on this situation. 

The general one-dimensional solution of eq. ([!]) is a parabola characterized by 
three parameters. In addition to the two diffusion equations, four boundary 
conditions are necessary to determine the two distributions p\ and p 2 ■ 



Pl {-£/2) = (2) 

Pl (£/2-R inc ) = p 2 (£/2-R inc ) (3) 

p[(£/2-R inc ) = p' 2 (£/2-R inc ) (4) 

-Dp' 2 {£/2) = y P2 (£/2) (5) 

<-i 



Condition (g) is for the special case of perfect absorbing step edges. (|3|) and 
OD) are necessary to obtain a smooth density between region 1 and 2. The 
left hand side of @ is the particle current at the step edge. On the right 
hand side this is reformulated using the number of jump attempts Dp 2 (£/2) 
multiplied by the probability of overcoming the Ehrlich-Schwoebel barrier 
Eg- This probability is expressed as the inverse of a typical length £\ 



Figure 2: Arrangement of steps around the bottom terrace of width £ 
between two mounds. 



where a stands for the lattice constant. 

The resulting density distribution has the form as indicated in fig. |l|: a 
parabola in the upper region and linear close to the downward step. The 
detailed expressions of p\ and P2 are not of much interest since the evolution 
of the crystal is determined by the currents at the edges. In the following we 
will call u(£) the upward current, i.e. 

u{£) = -Dp\{-£/2) 

= - 2a J + h) [l 2 + m - 2R inc h - R inc 2 ) . (7) 

The downward current due to diffusion (the contribution of the incorporation 
mechanism is not included) becomes 

d{£) = -Dp> 2 (+£/2) 



F 



{£ - Rn 



,2 



2a(£ + £ 1 ) y " ~ mcJ 

Note, that these results are very similar to the corresponding equations (2.2) 
and (2.3) of reference j|] where no incorporation was considered. Setting 
-Rinc = we regain their results. 

The absence of a dependence on D reflects the ansatz of a quasi-stationary 
distribution. All arriving particles are compensated for by the loss of particles 
at the borders and hence the currents are proportional to F . The density 
itself is proportional to the ratio F/D which again is intuitively clear. 

3 Closure of bottom terraces 

In the following, we will reinvestigate the discussion of concerning the 
closure of a bottom terrace (c.f. fig. 0). In the limiting case of an infinite 
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Ehrlich-Schwoebel barrier the dynamics of the steps become very simple. 
We denote by x(t) the position of the right step which of course will depend 
on the time t. The origin is chosen to be in the middle of the bottom 
terrace. Due to the infinite Ehrlich-Schwoebel barrier the movement of the 
right and the left step are symmetric. The evolution is then described by 
x(t) = —Fx(t) —FRi nc . The first term corresponds to the particles which do 
fall on the bottom terrace and diffuse to the right. The second term is the 
contribution of particles which are incorporated from the step above (which 
is valid as long as the bottom terrace is more than a distance Ri nc away from 
a top terrace). As a result x(t) evolves as 

x{t) = (x + R inc ) exp {-Ft) - R inc . (9) 

As long as R inc > there exists a closure time 

tc= i ln( xo + R^y 



F \ Rl_ 

Without an incorporation mechanism (R inc =0) the bottom terrace will never 
be closed. This is the reason why Elkinani and Villain called their model 
the Zeno-model to remind the greek philosopher and his paradox. Even 
though the situation is changed if the discrete structure of the terraces is 
considered^ they showed that this trend still holds which gives rise to the 
formation of deep cracks. Likewise they found that even finite values of the 
Ehrlich-Schwoebel barrier do not change this growth scenario which has been 



investigated in more detail in ||25|| . Once mounds are built up they remain 
forever with a fixed lateral size. Our discussion of this limiting case shows 
that the inclusion of an incorporation mechanism changes the growth in a 
fundamental manner. 



4 Growth dynamics 

To set up the basic ideas of the behaviour during crystal growth we show 
two typical surface profiles according to the numerical integration of the 
step system. In fig. ^ we compare the resulting structure of the Zeno model 
H] without an incorporation mechanism and with the inclusion of such a 
mechanism. 

1 The currents can be translated into probabilities of placing a particle at the step edge. 
Hence, a bottom terrace of width one always has a nonvanishing probability to be filled. 



6 



Figure 3: In (a) we show the typical height profile of the Zeno model. With 
inclusion of an incorporation mechanism in (b) we obtain structures with 
slope selection and prevent the formation of deep cracks. In both cases we 
have deposited (from bottom to top) 20, 50, 100, and 200 ML where the 
height profiles have been shifted by +25, 0, -45, and -140 ML respectively. 
The regions are subsections of 200 a out of a surface of size 485 a. 

The simulations were carried out on on a system of 485 a width with param- 
eters corresponding to the model used in sec. |6|. 

19 / 0.9 eV \ a 2 a 2 
= eXP ( + 5^) a=S8 ' 2a 

-Rinc = la 

F = 1 ML s" 1 

As in [|J the Ehrlich-Schwoebel barrier has been suppressed for bottom ter- 
races of one lattice constant width. Without an additional incorporation 



mechanism the appearance of trenches is unavoidable in accordance to p5 
The incorporation mechanism gives rise to a well defined slope which does 
not change with time. Another fundamental difference is the coarsening be- 
haviour. Without an incorporation mechanism the trenches are stable and 
the number of mounds remains constant. The additional incorporation mech- 
anism leads to a coarsening behaviour. 

In lattice models as well as for continuum equations the coarsening is driven 



by fluctuations |26], |24J] and in 1+1 dimensions the corresponding exponent 



is 1/3. This is in accordance to Ostwald-ripening which has been predicted 



from the similarities of the relevant continuum equations fll7j . However, 
since we treat the step evolution in a deterministic manner we do not obtain 
a scaling behaviour. The only way fluctuations come into play during the 
simulation is when new islands are nucleated. As a consequence the evolution 
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of e.g. the width of the height distribution w is characterized by jumps (data 
not shown). A jump in w appears each time when two mounds merge. These 
findings are a direct confirmation of the relevance of the fluctuations for the 
coarsening behaviour. 



5 Slope selection 

The inclusion of an incorporation mechanism leads to slope selection which 
is apparent from fig. |3|b. Siegert and Plischke |16|] required a cancelation of 



upward and downward currents in the continuum equations. Again, in the 
case of an infinite Ehrlich-Schwoebel barrier the calculations are straight- 
forward. In this case the downward current on a vicinal terrace of size £ is 
solely due to the incorporation mechanism, i.e. proportional to FR inc . All 
the remaining diffusing adatoms will contribute to the upward current and 
hence the current will be F(£ — -R mc ). As a consequence the slope selection 
will be achieved with a mean terrace width of size 

t = 2R mc (11) 



in accordance to the findings in [24] 



It remains to calculate the terrace widths for arbitrary parameters. Since we 
know the currents (equations (|7|), (|8]) and the incorporation mechanism) we 
obtain the overall slope (resp. terrace width) dependent current 

J(£) = u(£) + d(£) + FR inc 

= 2 ^ F +li) ( 2 ^mc 2 + 4flincll ~ 2l£t) (12) 

Note that a positive J(£) signifies a downward current (to the right in fig. [[]). 
The stable slope, where no net upward or downward current remains, is given 
by the condition J(£)\ e=e * = and yields 



2-Rinc ~l" 



D 2 



'i 

2 



2i? inc + ^ e - E s/k B T (13) 
a 



As can be seen from expression (0) the current is positive for small values 
of i and becomes negative for £ > £*. Hence £* is stabilized by the current. 
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The stable slope does not depend on the diffusion constant. However, it 
should be clear from the derivation, that in order to achieve slope selection 
the typical diffusion length should be much larger than I*. Otherwise the 
vicinal terraces would not proceed via step flow. Rather new nucleation 
events on the terraces would lead to a rugged surface structure. 



6 Solid— On— Solid model 

In order to verify the predicted importance of the incorporation mechanism 
we use computer simulations of the Solid-On-Solid (SOS) model on a simple 
cubic lattice. All processes on the surface are (Arrhenius-) activated processes 
which are described by one common prefactor u and an activation energy 
which is parameterized as follows: Eb is the barrier for surface diffusion, at 
step edges an Ehrlich-Schwoebel barrier E$ is added. However, this barrier is 



not added for a particle which sits on top of a single particle or a row |27|, |23 
Each next neighbour contributes E^ to the activation energy. Within this 
framework the diffusion constant becomes D = u ■ exp (—EB/ksT). 

Here, we concentrate on a particular set of parameters even though other 
parameter sets were used as well. We choose u = 10 12 s _1 , E B = 0.9eV, 
E^ — 0.25eV, and Eg = O.leV. This model was already investigated in 
p8| and reproduces some kinetic features of CdTe(OOl). The deposition of 
particles occurs with a rate F. The two simulations shown in fig. || are carried 
out on a 300 x 300 lattice at 560 K and started on a singular (flat) surface. 

In fig. H the resulting surfaces with and without the inclusion of the incor- 
poration mechanism are shown. Without an incorporation mechanism no 
slope selection occurs. Clearly, without incorporation the configuration of 
the towers remains unchanged whereas the inclusion leads to coarsening. 
The number of mounds diminishes with time. Hence, without an incorpora- 
tion mechanism no coarsening can be identified. We want to mention that 
it seems that at higher temperatures the attachment/detachment kinetics of 
atoms at step edges yields a coarsening effect (data not shown). However, 
still no slope selection has been observed. 

At first glance our findings contradict previous results obtained with a very 
similar model. Smilauer and Vvedensky obtained a formation of mounds with 
slope selection irrespective of the inclusion or exclusion of an incorporation 
mechanism jTJJ]. However, they implemented the Ehrlich-Schwoebel barrier 
in a different way. Rather than to hinder the jump over a step edge they im- 
pede the jump towards a step edge. Their motivation for this implementation 
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Figure 4: We compare the morphology of the surfaces without (left) and 
with (right) an incorporation mechanism. Note that the two grey-scales are 
different. The heights in the left picture range from 1208 to 1326 whereas 
the right surface only spans a height difference of 12 from minimum to max- 
imum. The contour lines are drawn for the same surface at an earlier stage 
where only 300 ML have been deposited. Without incorporation the mounds 
(towers) are nearly unchanged despite the deposition of 1000 ML. 



was to allow the adatoms to leave a small line of particles of width one which 



has been tested as a cause for reentrant layer- by-layer growth p7 , [23[| . In our 
simulations the same goal is achieved by suppressing the Ehrlich-Schwoebel 
barrier in such a situation. However, in their simulations particles arriving 
directly at a step edge have a probability of 1/4 to jump down the edge, 1/4 
to jump away from the edge and 1/2 to jump along the step edge. Effectively 
this leads to an incorporation radius of length 1/2. 



Other simulations of SOS-models used bcc(001) |29|, [15|] in order to study 
the growth of typical metals. In these simulations the SOS-restriction is 
implemented in such a way that an adatom must be supported by the four 
underlying atoms. Hence, the downward funneling process is directly im- 
plemented. Again, as a result slope selection is achieved, which has already 
been discussed in great detail in pO . 



7 Detachment and desorption 

One might assume that other mechanisms could lead to a zero in the slope 
dependent current. In the following we will carry out an analogous calcula- 



te 



tion with an adatom-detachment rate from steps and inclusion of desorption 



pl |. One might assume that both processes generate an effective downward 
current which can compensate for the Ehrlich-Schwoebel effect. To inves- 
tigate whether they are sufficient to obtain slope selection (and to simplify 
notation) we exclude the incorporation mechanism. Thus, the distinction of 
the two regions on a terrace is not necessary. 

The desorption of diffusing adatoms is easily incorporated including a term 
—p(x)/r in the diffusion equation ([!]) In order to include detachment 

from steps we have to replace boundary condition (^j by 

-D P \-e/2) = 1 --p(-e/2) (14) 

a 

where 7 stands for the detachment rate from steps. Accordingly, the bound- 
ary condition at downward step has to be corrected and reads now 

- Dp'{l/2) = £p(*/2) - 7^. (15) 

The overall slope dependent current becomes 

(A - 1) (it - a) - DF) 
m = ; JW7-. J U ~ } — - (16) 



(A + a) V^(A + 1) + ^(A — 1) + D(A — 1 

where 

A = e 2e ^ 

has been introduced. Note that A is always greater than one. 

To discuss the qualitative behaviour it is sufficient to look at the numerator 
of J(£) (the denominator is always positive). The first important result is 
that no slope selection is possible. Only for I = the current is zero (of 
course, there is no upward and downward current as well). 

Even though there is no slope selection, one can discuss whether growth will 
proceed via layer-by-layer growth (J(£) > for all t, i.e. terraces tend to 
grow larger) or a growth instability is present (J(£) < for all £, i.e. particles 
are preferably incorporated at the upper steps). 



In the well known limit of negligible detachment or desorption rates |32 
(7 — > or r — ► 00) the Ehrlich-Schwoebel effect alone determines the sign of 
J(£). As expected, for positive step edge barriers {t\ > a) growth becomes 
instable whereas negative values of Es stabilize layer-by- layer growth. 
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If we assume a positive E$ in the general case we obtain a critical current 

F„ = £ (17) 
where the current changes its sign. 

If one expresses the diverse rates as used in the Solid-on-Solid simulations 
and sets a = 1 one obtains 

F c = Ud e- {ED+E ^d)/kBT_ (lg) 

where the desorption rate u d e~ ED ^ BT has been introduced. E hind represents 
the typical binding energy of a detaching adatom. In 2+1 dimensions this 
should be approximately -Ebind ~ 2E N . Using the model of the previous 
section, v d = u Q) and E D = l.leV (parameters which are a reasonable guess 
for CdTe(OOl), J28[) one obtains a critical flux F c = 0.004 ML/s. However, 
this crossover should not be observable in experiments, since at such low 
external fluxes the step flow of the preexisting steps will dominate the surface 
evolution. 



8 Saturation profile with infinite step edge 
barrier 

In this section we will calculate the saturation profile for the model with 
infinite step edge barrier. The discussion of the closure of the bottom terrace 
already showed that in this limit the calculations become very simple. As 
for the bottom terrace the dynamics of higher steps become independent of 
the above lying terrace. The steps Xi(t) evolve according to 

dxi 



(t) = -F( Xl (t) + R mc ) 
(t) = -F(x 2 (t) - xx(t)) 



dt 
dx 2 

H 

fi t ■ 

-£(t) = -F( Xi (t) - xUt)) (19) 



Measuring the time in units of 1/F (i. e. setting F = 1 in the above equations) 
the time to grow one monolayer is equal to one. In addition, to simplify 
notation we will measure all lengths in units of R inc . 
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i 


a*(0) 


Xi (0) - x^O) 


1 


e- 1 


1.71828 


2 


e 2 - e - 1 


1.95249 


3 


e 3 - 2e 2 + |e - 1 


1.99493 


4 


e 4 - 3e 3 + 2e 2 - |e - 1 


2.00009 


5 


e 5 - 4e 4 + |e 3 - f e 2 + ^e - 1 


2.00007 



Table 1: Analytical expressions of the saturation profile step positions and 
the numerical values of the terrace widths. 

The solution for the lowest terrace has been given in eq. (H). The solution 
for all steps is 

n ^n— i 

x n(t) = E 7-^j( 1 + " 1 (20) 

i=l \ n % )- 

as can be easily verified. If we want to calculate the steady state saturation 
profile we have to require 

x i+1 (l) = Xi (0). (21) 

It should be stressed that this is the only assumption: the surface morphology 
is a self-reproducing structure. If the upper steps x i+ \ after deposition of one 
monolayer would be greater than Xj(0) this would result in a flattening of 
the surface. Otherwise the slope would become steeper. Using the solution 
for the bottom terrace © we obtain the initial value 

ar 1 (0) = e-l (22) 

when we require that the bottom terrace will be closed at time t = 1. 
For the upper terraces eq. (|20| ) yields a recursion relation 

n 1 + r .(Q) 

x n {\) + 1 = x n _!(0) + 1 = E 7 I, ( 23 ) 

which can be solved as described in the appendix using the generating func- 
tion. As a result one obtains the initial positions of the steps on an infinite 
symmetric step profile. Every time the bottom terrace is closed the steps 
(with new indices) are located at these positions. 

In table we show the analytical expressions for the step positions derived 
from the generating function. In addition, the numerical values of the ter- 
race widths are shown. With growing index the terrace widths are rapidly 
approaching two. Even though they oscillate around this value it can be 
shown that 

lim(a; i (0)-x i _i(0)) = 2. (24) 
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Note that we measure the lengths in units of i?i nc . Hence we predict a slope 
selection with slope l/(2i? inc ) as derived from the simple argument of section 

The derivation shows that the selected slope is controlled by the closure of 
the bottom terrace only. Another length scale is the nucleation length. It is 
defined by the typical length of a top terrace at which nucleation of an island 
occurs []33] . This length scale is responsible for the rounding of the towers in 
fig. ||b and causes a perturbation of the steady-state saturation profile. 



9 Conclusion 



We have investigated the effect of an incorporation mechanism on the mor- 
phology of growing surfaces. The inclusion of an incorporation mechanism 
in a 1+1 dimensional BCF-theory as well as in SOS computer simulation in 
2+1 dimensions is necessary in order to obtain slope selection and a coarsen- 
ing process. We were able to derive analytically the temperature dependence 
of the selected slope. We found that the Ehrlich-Schwoebel barrier alone 
controls the temperature dependence. In the limit of an infinite step edge 
barrier we derived the steady state saturation profile. In this case the result- 
ing mound morphology is controlled by the closure of the bottom terrace. 



A Generating function 

To simplify notation we introduce the shifted step positions bj = Xj(0) + 1. 
We will try to extract the generating function 

oo 

/(*) = 5>** (25) 

3=0 

for the shifted step positions. Clearly, the bj are only of physical meaning if 
j > and bo can be chosen arbitrarily. Starting from equation (|23|) 



n b- 

b n _ x = J2 7 l -TT- for all n > 2 (26) 



i=i (n-ty.e 
fri {n-i)\ 

oo oo n i n 

ez J2b m z m = EEt^aT (28) 



m=l n=2 ; 
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Choosing b = and using 61 = e we arrive at 



00 n t n 

e z f(z) = ' — ~ez 



n — i) 



n=0 i=0 

e -z f{z) = f(z)e z -ez 



f(z) 



(29) 
(30) 

(31) 

- ~ jl ^ (32) 

derived from the generating function are shown in table (§). In addition, the 
generating function can be used to formally prove equation (|24|). 



Thus, we finally obtain 



The lowest coefficients 



e z ~ l — z 



1 &f 
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